On the criterion for Bose-Einstein condensation for particles in traps 
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We consider the criterion for Bose condensation for particles in a harmonic trap. For a fixed angular 
momentum, the lowest energy state for a cloud of bosons with attractive interactions is the ground 
state of the cloud with all the angular momentum in the center-of-mass motion, and the one-particle 
reduced density matrix generally does not have a single large eigenvalue, but a number of them, 
suggesting that the state is an example of a fragmented condensate (Wilkin, Gunn, and Smith, Phys. 
Rev. Lett. 80, 2265 (1998)). We show that a convenient way to describe correlations in the system 
is by defining an internal one-particle reduced density matrix, in which the center-of-mass motion 
is eliminated, and that this has a single eigenvalue equal to the number of particles for the problem 
considered here. Our considerations indicate that care is necessary in formulating a criterion for 
Bose-Einstein condensation. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40.Db, 67.40.Vs 
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In investigations of the properties of the helium liquids, 
rotation has proved to be a valuable probe, and this has 
led to the discovery of vortex lines in liquid He and in 
the superfluid phases of liquid 3 He. In the discussion of 
superfluidity and Bose-Einstein condensation it is usu- 
ally assumed that particles are confined in a box, with 
a uniform potential throughout its interior, and with in- 
finitely high walls. This situation is very different from 
that encountered in studies of dilute atomic vapors in 
traps, where the potential is generally closer to a har- 
monic oscillator one. We therefore consider the proper- 
ties of rotating atomic clouds of bosons at zero temper- 
ature. Among possible ways of putting angular momen- 
tum into the system are the creation of vortex lines 
and excitation of surface modes and other internal exci- 
tations g . For traps with an axis of symmetry, another 
possibility is to excite the center-of-mass motion, which 
for modes having rotational quantum number m z = ±1 
about the symmetry axis is equivalent to excitation of 
the simplest surface mode. The purpose of this paper is 
to consider the criterion for Bose condensation in a cloud 
of particles in a trap. 

A criterion for Bose-Einstein condensation was pro- 
posed by Penrose ||, and subsequently elaborated by 
Penrose and Onsager 0], and by Yang JtJ. Consider the 
one-particle reduced density matrix. This is defined in 
terms of the many-body wave function ip(ri, ...r^v), where 
the Vi are the coordinates of the particles, by the equation 

p (1) (r, r') = J dr 2 . ..dr N ip(r,r 2 ■ . . r N )yj* (r', r 2 . ..r N ). 

(1) 

The density matrix may be expanded in terms of its 
eigenfunctions Xji r ) with eigenvalues Xj\ 



pW(r,r')=^A jX *(r')x J -(r). 



(2) 



The condition generally adopted for the existence of Bose 
condensation is that there should be one eigenvalue that 
is of order the number of particles, in the limit when the 
number of particles tends to infinity. However, as we 
shall indicate below, this can be misleading. 

In an instructive paper Wilkin et al. [pi considered 
bosons in a harmonic trap, and showed that for weak 
attractive forces the lowest state with a given angular 
momentum is one in which all the angular momentum 
is associated with the center-of-mass motion. This cor- 
responds to excitation of a dipolar (magnetic quantum 
number m z = ±1) surface mode of oscillation ||. An 
analysis of the eigenvalues of the one-particle density ma- 
trix revealed that there is not a single large eigenvalue, 
but rather there can be many eigenvalues of a compara- 
ble size. This led the authors to suggest that the state of 
the rotating system corresponds to that of a fragmented 
condensate of the type proposed by Nozieres and Saint 
James |jTo|| . However, this conclusion seems surprising, 
because the state is just the ground state of the cloud, 
which is Bose-Einstein condensed, executing center-of- 
mass motion. The problem is to identify a quantity that 
can reveal the Bose condensation. 

Let us consider N spinless bosons in a harmonic trap. 
This system has the important property that the center- 
of-mass and relative motions of the particles are separa- 
ble, and therefore the wave function may be written in 
the form 



r/>(ri, ...r N ) = ^cm(R-)V'rel(qi,q2, ■ ■ ■ ,qjv)> 



(3) 



where the center-of-mass coordinate is given by R = 
(l/iV)E^ =1 rj and = r - ; — R is the coordinate of particle 
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i relative to the center of mass. Note that the definition 
of the center-of-mass coordinate implies that only N — 1 
of the are independent. Let us now introduce a quan- 
tity which we call the internal density matrix 

Pintfoi, = / dR^dcu + R})V*(K + R}) 

= Vrel({qi})d({qa)- (4) 

This has the property that the center-of-mass motion has 
been eliminated, and consequently it depends only on co- 
ordinates relative to the center of mass. One may also 
construct reduced density matrices in terms of the rela- 
tive coordinates in the standard way. For example, the 
internal one-particle reduced density matrix is 

Pint(qi>q'i) = J rfq2"-rfqjv-iPint(qi,q'i)- (5) 

Consequently the internal density matrix for the lowest 
energy state of a cloud of bosons with weak attractive in- 
teractions and nonzero angular momentum in a harmonic 
trap is the same as for the ground state in the absence of 
rotation. It has one large eigenvalue which is equal to N, 
to within terms of unity, the deviation from N being due 
to the constraint that the sum of the coordinates relative 
to the center of mass must vanish. 

Our discussion indicates that caution is required in for- 
mulating a criterion for Bosc-Einstcin condensation in 
finite systems, and that the one generally used can be 
misleading. For systems in harmonic traps the internal 
density matrix introduced above can be a useful concept, 
because it removes the effects of the center-of-mass mo- 
tion which, for a rotating cloud, tends to spread out the 
weight of the usual one-particle reduced density matrix 
over a number of states. 

We are grateful to Ben Mottelson for helpful conversa- 
tions. 



momentum the ground state for a given angular momen- 
tum is obtained for weak coupling (Na/a osc <C 1, where 
a is the scattering length for two-body collisions and a osc 
is the usual characteristic length scale for the harmonic 
oscillator) by adding quadrupolar (m z — ±2) surface ex- 
citations to the non-rotating ground state (B. Mottel- 
son, Phys. Rev. Lett. 83, 2695 (1999), G. M. Kavoulakis, 
B. Mottelson, and C. J. Pethick, (in preparation)). As 
the repulsive coupling increases, the optimal excitations 
become surface modes with higher and higher values of 
\m z \. The nature of the lowest energy state for a fixed 
non-zero angular momentum therefore changes discon- 
tinuously as the coupling passes through zero. 
[10] P. Nozieres and D. Saint James, J. Phys. (Paris) 43, 1133 
(1982). 
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